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Abstract 

A fundamental issue that limits the efficiency of many photoelectrochemical systems is 
that the photon absorption length is typically much longer than the electron diffusion length. 
Various photon management schemes have been developed to enhance light absorption; one 
simple approach is to use randomly scattering media to enable broadband and wide-angle 
enhancement. However, such systems are often opaque, making it difficult to probe photo- 
induced processes. Here we use wave interference effects to modify the spatial distribution 
of light inside a highly-scattering dye-sensitized solar cell to control photon absorption in a 
space-dependent manner. By shaping the incident wavefront of a laser beam, we enhance or 
suppress photocurrent by increasing or decreasing light concentration on the front side of the 
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mesoporous photoanode where the collection efficiency of photoelectrons is maximal. En¬ 
hanced light absorption is achieved by reducing reflection through the open boundary of the 
photoanode via destructive interference, leading to a factor of two increase in photocurrent. 
This approach opens the door to probing and manipulating photoelectrochemical processes in 
specific regions inside nominally opaque media. 
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Introduction 

Optieal seattering, indueed by struetural disorder, ean inerease the dwell time of light inside an 
absorbing material, thus enhaneing the ehanee of absorption.-”- Compared to other light-trapping 
sehemes that have been developed for solar eells, the non-resonant nature of this seheme enables 
broadband enhaneement for light with a wide angle of inoidenoe.->^ However, the degree of absorp¬ 
tion does not depend solely on the absorption eoeffieient and seattering properties of the medium. 
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but also on the coherence and spatial characteristics of the input light. One recent example con¬ 
firming this is the demonstration of coherent perfect absorption, which is the time-reversed analog 
of laser emission.-^ Similar to lasing, this effect is narrow-band and requires tuning of both the 
frequency of the input light and absorptivity of the material to achieve 100% absorption.—-— A 
related but distinct approach, termed coherently enhanced absorption, takes advantage of the large 
number of degrees of freedom available for multiple scattering in a lossy environment to control 
the absorption.-i^*i^ It has been shown theoretically an appropriate coherent superposition of in¬ 
put fields can significantly enhance the coupling of light into a scattering medium, increasing the 
absorption at any frequency over a broad range. 

Recent advances in wavefront shaping (WFS) techniques using spatial light modulators (SLMs) 
have enabled the generation of arbitrary wavefronts at the input of an optical system. This tech¬ 
nique is driving forward the development of optical imaging techniques through scattering media 
such as biological tissue and fog,—-— and also enables the control of light intensitypo¬ 
larization state,— >2^ spectrum,— or temporal profile^-— at the scale of wavelength (single or a 
few speckles). While the wavelength-scale control is desired for imaging with high resolution,— 
the ability to manipulate light distributions on a much larger scale in a disordered system is needed 
for many applications. One example is photoelectrochemical cells, which often utilize light scat¬ 
tering to increase photon dwell time to enhance absorption. As the photoelectrons are generated 
in a volume that contains hundreds of thousands of speckles, increasing the intensity of a single 
speckle even by a few thousand times makes little difference to the total photocurrent. Hence, 
global control of the light distribution is necessary, but this is known to be a much more challeng¬ 
ing optimization problem than local control. Good progress has been made recently to enhance the 
total transmission through a lossless opaque medium via WFS of an incident laser beam.— >^-— 
However, coherent control of light distribution and absorption over a macroscale region in highly 
scattering media has never been demonstrated. 

Here we use WFS to both enhance and suppress the photocurrent in a dye-sensitized solar 
cell (DSSC), with a photoanode comprised of a highly scattering layer of dye-sensitized Ti02 
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nanoparticles. Since the collection efficiency of charge carriers is highly dependent on where 
light is absorbed,— it is advantageous to trap light inside the scattering layer as close as possible 
to the interface with the electrode, where charges can be efficiently separated, thus minimizing 
charge recombination.— However, strong light leakage through the interface makes it extremely 
difficult to confine light to a region close to the open boundary and enhance the absorption there. 
But here we show a significant enhancement of photocurrent by optimizing the wavefront of an 
incident laser beam, and illustrate the underlying mechanism via numerical modeling. In addition, 
we show that WFS can also suppress absorption by increasing the amount of light reflected. This 
work demonstrates the great potential of WFS toward depositing energy into targeted regions inside 
turbid media, which will enable the spatially resolved characterization of light absorption and 
charge transport inside regenerative and artificial photosynthetic systems. 


Control of DSSC photocurrent 

A DSSC comprises a porous layer of sintered Ti02 nanoparticles covered by a molecular dye, in 
this case the ruthenium-based N719 dye,— deposited on a conductive glass substrate. It functions 
as the light-absorbing working electrode of the solar cell, while a platinum-loaded conductive 
glass slide acts as the counter electrode. An electrolyte fills the gap between the two electrodes 
and shuttles charges between them. The working electrode of a DSSC constructed in this manner 
can be considered as a scattering sample with homogeneous absorption.— Figure [Ik presents a 
schematic of the solar cells created in this study. Light is incident onto the porous Ti02 layer 
through the glass substrate coated with transparent conducting oxide. After light absorption, the 
excited dye molecules rapidly inject electrons into the Ti02. These electrons diffuse through the 
sintered particle network to be collected at the front side electrode, while the dye is regenerated 
via reduction by a redox mediator in the electrolyte, in this case the iodide/triiodide redox couple. 
The photo-induced current from a DSSC is typically measured without applying a bias voltage, 
and is known as the short-circuit current Isc- Not all photo-generated electrons can diffuse through 
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Figure 1: Experiment setup, (a) Sehematie of the dye-sensitized solar eell (DSSC) used in this 
study, whieh employs a thinner absorbing layer and a different eleetrolyte from the eonventional 
DSSC to demonstrate eoherent eontrol of photovoltaie eurrent. (b) Sehematie of the experimental 
setup (see Methods). A laser beam at A = 532nm is split to two orthogonal polarizations and 
projeeted onto the SLM. The phase front of both polarizations are modulated and then reeombined 
at a beam splitter. The shaped wavefront is imaged onto the surfaee of the porous eleetrode in the 
DSSC using a mieroseope objeetive. The illuminated area is a disk of diameter ~ 8.3 /im. The 
short-eireuit eurrent of the solar eell is measured using a potentiostat. 
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the Ti 02 layer and reaeh the eleetrode, beeause eleetrons may recombine or get trapped on the 
way.— The farther the electrons are generated from the electrode, the lower their probability of 
reaching it and contributing to the current. Hence, Isc depends not only on the amount of light 
being absorbed but also on the position where the light is absorbed inside the Ti 02 layer.— 

In the WFS experiment, the solar cell is excited by a narrow band laser beam at wavelength A 
= 532 nm. The input wavefront is modulated using the experimental setup schematically shown 
in Fig.[Tj5. The two linear polarizations are modulated independently— (see Methods) and recom¬ 
bined into a single beam, which is projected onto the scattering medium in the solar cell using a 
microscope objective (NA = 0.95). We measure Isc generated by the DSSC for each input phase 
mask and use a genetic algorithm to find the wavefronts that maximize or minimize Isc- The ge¬ 
netic algorithm is known to be robust against noise in WFS experiments.—!^ A cost function rfe, 
which estimates the variation of the conversion efficiency of the solar cell, is defined as: 

{he ~ I dark) / ^0 

where (Jsc) and (<^o) are the average short-circuit current and input optical power, respectively, 
measured for a set of random phase masks. Idark is the small dark current measured at = 0/iW 
{SI). 

We constructed the DSSC using standard procedures (see Methods and SI). The transport mean 
free path of light in the 20 ± 1 /im thick Ti 02 layer, I, is measured in a coherent backscattering 
experiment to be 3.2 /im.— To shorten the response time of the solar cell in order to perform WFS 
experiments on reasonable time scales, we use a concentrated non-volatile iodine-based electrolyte 
that also absorbs light at 532 nm. The absorption of light by the electrolyte is proportional to 
that by the dye at any location, as both scale linearly with local field intensity. Thus the short- 
circuit current Isc, which results from light absorption by the dye, also reflects the absorption of 
the electrolyte, which does not contribute directly to the current. The ballistic absorption length 
la is about 65 /im, and the diffusive absorption length (the typical distance light will travel via 
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Figure 2: Control of the DSSC eurrent and power by shaping the input wavefront of a laser beam, 
(a) The measured eost funetion r/g for the short-eireuit eurrent, normalized to the value with random 
input wavefront (green dashed line), inereases (blue solid line) or deereases (red dotted line) during 
the optimization of the input wavefront by a genetie algorithm. The photo-indueed eurrent is 
enhaneed 2.1 times or suppressed 3.0 times by wavefront shaping, (b) Output eleetrieal power 
of the DSSC illuminated by a random phase pattern (dashed green line), and the two optimized 
wavefronts for maximal rje (solid blue line) and minimal rje (dotted red line). The ineident laser 
power is fixed at 100/iW. The asterisk indieates the maximal power, whieh is enhaneed 4.7 times 
orredueed 5.7 times. 


multiple seattering prior to being absorbed). La = y/ffa/S, is 8.4 /im.—>^ Without absorption, the 
total transmission T of the Ti02 layer is estimated to be 26%, and the refleetion R = 74%. In the 
presenee of absorption, T is redueed to 7%, and R to 47%, with the remaining 46% being absorbed 
(in the absenee of wavefront optimization). These estimations agree well with the experimental 
measurement (SI). 

By optimizing the ineident wavefront, the absorption eost funetion rje is gradually inereased 
or deereased as shown in Fig. [2^. With the optimal wavefronts, the photoeurrent is enhaneed by 
a faetor of ~ 2.1 or suppressed by a faetor of ~ 3.0. We are thus able to eoherently vary the 
energy eonversion effieieney by a ratio of ~ 6.3. To quantify the variation of the output power 
due to the modifieation of the input optieal pattern, we measure the eurrent-voltage (I-V) eurves 
for the two optimized wavefronts eorresponding to maximum and minimum rjg, as well as for a 
random illumination pattern. Figure Wp displays the light-indueed eleetrieal output power, given 
by (I — hark) y, as a funetion of bias voltage for eaeh seenario. The maximum output power is 
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enhanced 4.7 times or reduced 5.7 times, leading to total variation by a factor of 27. 

Since approximately 46% of the light is absorbed in the solar cell for a random illumination, the 
highest possible enhancement of the absorption would be a factor of 2.2. If I sc were proportional 
to the total amount of light being absorbed,— the experimental enhancement of 2.1 of the short 
circuit current would correspond to an almost perfect absorption of light in the absorbing medium. 
However, achieving the maximum absorption requires complete control over the amplitude and 
phase of the light in all input channels.— Experimentally we control only a limited fraction of 
the input channels due to the finite numerical aperture of the objective lens and the fact that the 
probe beam is incident only on one side of the scattering sample. Moreover, our SLM setup only 
modulates the phase and not the amplitude of the incident light, as would be necessary to achieve 
the global optimum. Consequently, the enhancement of total absorption that could be achieved 
with WES is significantly reduced.— 


Manipulation of local absorption 

To investigate the cause for such a high enhancement factor measured experimentally, we perform 
numerical modeling. Without loss of generality, we simulate the scattering medium confined in a 
2D waveguide with reflecting sidewalls to reduce the computation load; the system dimensions, 
as well as the amount of scattering and uniform absorption, are chosen to be near to the experi¬ 
mental values (see Methods). Eight with transverse-magnetic polarization {E field perpendicular 
to the waveguide plane) is injected from one end of the waveguide, and its propagation in the ran¬ 
dom system is calculated with the recursive Green’s function method.— The transmission and 
reflection of light for the systems with and without absorption, averaged over random input wave- 
fronts, match the measured values of real samples. To model the effect of WES in our experiment, 
we use the empty waveguide modes as the basis to shape the input field. During the optimization 
process, we modulate the phase of each empty waveguide modes while keeping their amplitude 
constant. The (linear) space-dependent absorption of light is proportional to the local intensity 
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\E{x^y)\'^, and the total absorption is proportional to the spatially integrated intensity. To take into 
aeeount the spatial variation of the eleetron eolleetion effieieney D{x) = exp[—where ^ is the 
eleetron diffusion length, we define a new eost funetion rj^ whieh models the relative ehange of the 
short-eireuit eurrent similarly to r/g: 

11 Dix)\E{x,y)\^dxdy 
il lD{x)\Eo{x,y)\^dxdy)' 

where E{x^y) is the eleetrie field distribution in the seattering medium for the optimized input 
wavefront, EQ{x,y) is the field distribution for a random input wavefront, and (..) represents the 
average over random input wavefronts. As in the solar eell, the absorption of light that oeeurs 
far from the front surfaee of the random medium eontributes little to rjs- We used both a genetie 
algorithm and a sequential algorith m to maximize or minimize rjs, and found similar results, 
although the sequential algorithm eonverges faster than the genetie algorithm due to the absenee 
of noise in simulation. 

Figure [3^ shows the enhaneement of the eost funetion rjs in the random medium with a uniform 
absorption eoeffieient. For eomparison we also simulate the ease of D{x) = 1, whieh eorresponds 
to a uniform eost funetion similar to the standard eoherent enhaneement of global absorption.-^^*^ 
Both the non-uniform {D{x) = exp[—x/(^]) and uniform eost funetions inerease eontinuously in the 
optimization proeess and eventually saturate; however as expeeted, the final enhaneement faetor 
for D{x) = exp[—x/(^] signifieantly exeeeds that for D{x) = 1. Figure[3b represents the integrated 
intensity I{x) = f \E{x,y)'\^dy inside the random medium. With random wavefront illumination, 
I{x) deeays exponentially in the disordered system, as expeeted for light diffusion in the presenee 
of absorption.— When the input wavefront is optimized for uniform eost funetion {D{x) = 1), I{x) 
is peaked in the interior of the random system. The buildup of energy deep inside the seattering 
medium enhanees the total absorption from 46% to 82%. Meanwhile, the refleetion deereases from 
47% to 9%, and the transmission slightly inereases from 7% to 9%. This result eonfirms that the 
global absorption is enhaneed by WFS through eausing the light to penetrate deeper into the lossy 
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Figure 3: Numerical simulation of WFS experiment, (a) The cost function rjs, defined in Eq.(2), 
increases and then saturates when the input wavefront is optimized for D(x) = exp[—x/i^] (red 
solid line) and D(x) = 1 (blue dashed line) using a sequential algorithm. ^ is taken to be 3 /im 
(see Methods), (b) The cross-section integrated intensity in the random system illuminated by 
random input wavefronts (black dotted line) and the optimized wavefront obtained with the non- 
uniform cost function for D{x) = exp[—x/(^] (red solid line) and the uniform cost function for 
D{x) = 1 (blue dashed line). The light is concentrated near the front side of the random system 
for D{x) = exp[—x/(^], but penetrates deeper into the system for D{x) = 1. (c,d) Same as (a,b) 
for the minimization of For D{x) = exp[—x/(^], the light intensity near the front side of the 
random system is reduced, but increased deeper into the system. For D{x) = 1, the light intensity 
is reduced in the interior of the system but remains almost the same near the front surface. 
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disordered medium. 


In contrast, the optimization for the non-uniform cost function {D{x) = exp[—x/i^]) leads to a 
significant increase of light intensity close to the front surface of the random medium [Fig. [3b], 
since only this region contributes significantly to the cost function rjs- Naively, the proximity to 
the open boundary would be expected to facilitate light escape from the absorbing medium, thus 
reducing the absorption and increasing the reflection. Instead, for the optimal wavefront, the total 
absorption is increased from 46% to 67%, and the reflection decreased from 47% to 28%. This 
result suggests some partial trapping of the light near the surface, leading to a buildup of its in¬ 
tensity in the vicinity of the open boundary where the front side electrode is located. Hence, the 
numerical simulation illustrates that our experimental procedure, which does feedback optimiza¬ 
tion based on the short-circuit current, does not simply optimize total absorption, but optimizes 
absorption locally, near the working electrode, so as to preferentially enhance the collected cur¬ 
rent. The measured current increase by a factor of 2.1 does not imply perfect absorption of all 
incident light and can be realized with the incomplete experimental control. 

Modal interference effects 

To understand how the leakage rate of the light close to the open boundary is reduced, we inves¬ 
tigate the contribution of different quasi-modes in the system to the reflected light. It has been 
shown previously that light propagation in random media can be well explained in terms of the 
interference effects of quasi-modes.— We take this approach by decomposing the electric field 
inside the random medium in terms of the quasi-modes of the system: 

E {x, y) = ainUm (v, y) + 52 ^mVm {x, y), (3) 

m m 

where Um (Vm) is the m-th mode with purely outgoing (incoming) wave boundary condition, and 
Urn = V*,.— The reflected field at the front surface {x = 0) of the disordered medium can be ex¬ 
pressed as Er{y) = Lm^mMm(0,y)- The total intensity of reflected waves is R = j\Er{y)\^dy = 
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where = amjE*{y)um{0,y)dy represents the eontribution of the m-th mode to the 
refleetion. While is a eomplex number, 7? is a real number, and the eontributions of dif¬ 
ferent modes to R interfere with eaeh other. The interferenee effeet ean be quantified by C = 
I Em C(m P/ Em I P (SI)- If C > 1 (C < 1), the modes interfere eonstruetively (destruetively) at the 
front surfaee to enhanee (suppress) refleetion. 
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Figure 4: Modal analysis of the wavefront optimization. Contributions a,„ of individual modes m 
supported by the random system to the total intensity of refleeted waves. The modes are exeited 
by a random input wavefront (a), the optimized wavefront for maximal rjs (b) and minimal rjs (e) 
with D(x) = exp[—x/(^]. In (a) the modes have random phases and their interferenee is averaged 
out, C = 1.05. In (b) the modes interfere destruetively, C = 0.21 < 1, to reduee refleetion and 
enhanee absorption. In (e) the modes interfere eonstruetively, C = 1.47 > 1, to enhanee refleetion 
and reduee absorption. 


To reduee the eomputing time for the full modal analysis, the length and width of the system 
in the simulation is redueed by a faetor of two. The seattering and absorption parameters are 
adjusted so that the total transmission, refleetion and absorption for random input wavefronts are 
essentially unehanged. We use the finite-differenee frequeney-domain (FDFD) method to ealeulate 
the quasi-modes whose frequeneies are elose to the frequeney of the input light. The ealeulated 
values of for the different seenarios with D(x) = exp[—x/(^] are plotted in the eomplex plane 
in Figs. |4h-e. The modes have uneorrelated phases for random wavefront illumination and their 
interferenee is averaged out (and therefore C is elose to 1). We see that the optimized wavefront 
for enhaneed loeal absorption eauses a destruetive interferenee among the modes, making C mueh 
smaller than 1, and more importantly, that refleetion is minimized and absorption is maximized at 
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the front surface. Despite the complexity of the medium, WFS exploits intereferences to efficiently 
concentrate light in the targeted region of the system. 

Similarly, we show that WFS can also reduce the energy in a chosen region by minimizing the 
cost function rjs [Figs. [St and d]. For D{x) = 1, the reflection is greatly enhanced, decreasing 
the light intensity deep inside the system and thus the total amount of absorption. For D{x) = 
exp[—x/<^], the minimization of the cost function is less efficient as it requires us to simultaneously 
increase the reflection and decrease the absorption close to the front surface. By decomposing the 
fields into the quasi-modes of the system, we confirm that the increase in reflection is caused by the 
constructive interference of excited modes [Fig. |4t], which facilitates light escape and suppresses 
absorption. 

Conclusion 

In summary, we have experimentally demonstrated coherent control of light distribution and ab¬ 
sorption over a length scale much larger than the wavelength of light inside an opaque mesoporous 
photoelectrochemical system that utilizes multiple scattering to improve light confinement. Our 
proof-of-concept experiment shows that it is possible to vary the photocurrent in a dye-sensitized 
solar cell (DSSC) by tailoring the input wavefront of a laser beam. To maximize the DSSC’s short- 
circuit current, the optimal wavefront concentrates the input light to the front side of the porous 
photoanode, where photo-generated electrons are most efficiently collected. Despite the proximity 
to an open boundary, we find that destructive interference among the excited modes minimizes 
light leakage and enhances absorption. 

Although in the current WFS experiment a narrow-band laser is used for coherent control 
of absorption, a recent theoretical study suggests that this approach is applicable to a relatively 
broadband light.— More generally, our method, assisted by a guidestar,— will allow for spatially- 
resolved characterization of photo-induced processes inside opaque media, such as photoelectro¬ 
chemical cells. This approach opens the possibility of investigating photoelectrochemical pro- 
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cesses in targeted regions inside dense seattering media, sueh as that used in regenerative and 
artifieial photosynthetie solar cells. Those studies will provide physieal insight into the operation 
of photovoltaie and photoeleetroohemieal systems, potentially leading to new methods to improve 
photoeonversion performanee. 

Methods 

Sample preparation. To faeilitate the WFS experiment, we made a DSSC with a thieker Ti02 
seattering layer, redueed loading of N719 dye, and a thin (175 /im thiek) tin-doped indium oxide 
(ITO) eoated glass substrate for the working eleetrode. The porous Ti02 working eleetrode is made 
from ~ 20nm Ti02 nanopartiele paste (P25, Degussa). The paste was doetor-bladed onto a 175 /im 
thiek ITO-eoated glass eoverslip (SPI, Produet No. 06477B-AB) 5 times. The sample was heated 
to 200°C for 10 min then eooled to room temperature between each application of a layer of paste. 
The final thiekness of the Ti02 layer is 20 ± 1 /im. The working eleetrode is sintered at 450°C for 1 
hour, then immersed in a solution of N719 (ditetrabutylammonium eis-bis(isothioeyanato)bis(2,2’- 
bipyridyl-4,4’-diearboxylato) ruthenium(II) for 30 minutes.— The eounter eleetrode is made by 
heating hexaehloroplatinie aeid in ethanol (2 mg/mL) to 400°C for 15 minutes on a fluorine-doped 
tin oxide (FTO) eoated glass slide. The solar eells are assembled with a 25 /im Surlyn spacer hot- 
pressed between the eleetrodes and filled with a non-volatile iodine-based eleetrolyte (Mosalyte 
TDE-250, Solaronix SA, Switzerland). 

Optical setup. A shown sehematieally in Fig. [H the output beam from a eontinuous-wave 
Nd:YAG laser is expanded (lenses L’ and L”) and projeeted onto a phase-only SLM (Hamamatsu, 
X10468-01). A half-wave-plate (Ai) is used to balanee the light intensities of two orthogonal linear 
polarizations. To eontrol independently the two polarizations, a polarizing beam splitter (PBS) 
attached to a right-angle prism (P) separates the two polarizations of the laser beam, whieh are 
modulated by different areas of the SLM. The eube and the prism are mounted in the same holder to 
eliminate fluetuations from independent refleetive elements in the two paths. A seeond half-wave- 
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plate (A 2 ) rotates the polarization to mateh that of the SLM. The two illuminated areas on the SLM 
corresponding to the two polarizations of the input light are divided into 1740 macropixels. The 
modulated wavefronts are then recombined and projected onto the pupil of a microscope objective 
(MO) with numerical aperture 0.95. A beam splitter (BS 2 ) along with a photo-detector (PD) are 
used to measure the optical power ,^0 impinging on the solar cell. The photo-induced current of 
the DSSC is measured by a potentiostat (Digi-Ivy DY2116b). 

Numerical simulation. In the 2D model, dielectric cylinders with radius 50 nm and refractive 
index 2.5 are randomly positioned inside a waveguide of width VP = 20 /rm and refractive index 
1.5. The length of the random structure is L = 20 jum, and the effective index of refraction is 
1.9. The wavelength of input light is A = 532 nm, and the number of transverse guided modes 
in the waveguide is 114. The transport mean free path £ is 3.7 jum. Without absorption and 
wavefront optimization, the transmission is T = 26% and the reflection is = 74%. A uniform 
absorption coefficient is introduced via the imaginary part of the refractive index to the random 
system. The diffusive absorption length is 8.5 jum. T is reduced to 7%, and reflection to 47%, thus 
the absorption A = 46%. These numbers match the measured values of our DSSC. The electron 
diffusion length ^ is chosen to be 3 /rm, close to the value of DSSCs similar to ours.— 

In the modal analysis to decompose the field distribution in the random medium with the quasi¬ 
modes, we keep the wavelength A the same, but reduce the number of waveguide modes to 37. £ 
is reduced to 2 jim, and the diffusive absorption length to 4.2 /im. Without absorption, T = 24% 
and R = 76%; with absorption, T = 1%, R = 46%, and A = 47%. The electron diffusion length ^ 
is also shortened to 1.5 /rm. 

Supporting Information Available 

More details on the synthesis and optical characterization of the DSSC. Additional numerical sim¬ 
ulation results of WFS to maximize and minimize the cost function rj^. Detailed description of 
modal analysis on the electric field distribution inside the random medium. This material is 
available free of charge via the Internet at http : / /pubs .acs.org/, 
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